A reconstruction from small-angle neutron scattering measurements of the real space 
magnetic field distribution in the mixed state of Sr 2 Ru04. 
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We have measured the diffracted neutron scattering intensities from the square magnetic flux lattice 
in the perovskite superconductor S^RuCU, which is thought to exhibit p-wave pairing with a two- 
component order parameter. The relative intensities of different flux lattice Bragg reflections over 
a wide range of field and temperature have been shown to be inconsistent with a single component 
Ginzburg-Landau theory but qualitatively agree with a two component p-wave Ginzburg-Landau 
theory. 



PACS numbers: 61.12 Ex, 74.60 Ge, 74.70 Tx 

The discovery of superconductivity at temperatures 
near IK in strontium ruthenate (lj has excited great in- 
terest because it is a superconducting layered perovskite 
which does not contain copper. However it shows great 
differences from the High-T c cuprates: it is a stoichio- 
metric undoped compound with a long mean free path, 
in which the electrons form a Fermi liquid with a well- 
established quasi- two-dimensional Fermi surface Q . Fur- 
thermore, it was suggested || that the strongly interact- 
ing electrons pair in a triplet p-wave state (rather than 
the singlet, mainly d-wave state which is believed to oc- 
cur in hole-doped cuprates). Clear evidence of non s-wave 
pairing in this compound has been provided by the ob- 
servation Q that nonmagnetic impurities strongly sup- 
press T c , which extrapolates to « 1.5 K in the clean limit. 
Strong support for triplet (p-wave) pairing is given by the 
results of Ishida et al. who have measured the Knight 
shift with a field parallel to the Ru02 planes ||; the spin 
susceptibility measured by the Knight shift is not sup- 
pressed below T c , unlike a singlet superconductor. Also, 
/iSR measurements in the Meissner state in zero field |(J 
have revealed spontaneous fields, which can be generated 
by domain boundaries, surfaces and impurities in a su- 
perconductor which breaks time-reversal symmetry 0. 
Such states can arise most naturally with p-wave pair- 
ing, but also are possible with d-wave singlet pairing. 

Agterberg |||9| argued that if the pairing was time- 
reversal symmetry breaking p-wave, then in tetragonal 
symmetry the d- vector ]l0| has the symmetry z-exp(±i<^) 
(cp is the azimuthal angle about the tetragonal c axis), 
and a two-component Ginzburg Landau (TCGL) the- 
ory would be expected to describe the superconductor. 
In zero field, this gives two degenerate states which are 
related by time reversal; with a field applied in the c- 
direction perpendicular to the planes, one is dominant, 
but the other is also present pH |. Under these conditions, 
a small amount of anisotropy in the Fermi surface would 
lead to a square flux lattice instead of a triangular one, 



with the orientation of the square flux line lattice (FLL) 
relative to the crystal axes determined by the orienta- 
tion of the fourfold anisotropy of the paired electrons. 
The FLL structure has been observed in this material 
Jl2[| and is observed to be square over a wide range of 
field and temperature. The nearest-neighbour directions 
in the square FLL are at 45° to the Ru-O-Ru directions 
in the crystal lattice |l3| . 

These results are consistent with the pairing wavefunc- 
tion described above. However, a square FLL is also seen 
in borocarbide superconductors, which are definitely non- 
p-wave ]l4|-|l6|]. Also, one can measure spontaneous fields 
in a superconductor by /iSR due to other causes or from 
other states than that proposed, and application of a 
strong field in the basal plane to observe the Knight shift 
might alter the pairing state. Hence, it is important to 
obtain further evidence as to what kind of superconduc- 
tivity occurs in strontium ruthenate. Here we present a 
detailed study of the scattered neutron intensities from 
the FLL. We show that they are not consistent with a sin- 
gle component Ginzburg-Landau model. Also we demon- 
strate how the local B(r) may be reconstructed from our 
data and show that the FLL structure is quite different 
from the Abrikosov one. 

We shall present measurements of intensities of higher- 
order Bragg reflections from the FLL so we consider how 
they are related to the FLL structure. The formula [l?]] 
for the integrated intensity Ihk of a (ft, k) diffracted peak 
of wave vector c\hk gives: 

hk oc (1) 
Qhk 

where Fhk is a spatial Fourier component of the local 
field B(y) in the mixed state: 

B(r) =^F^exp(iq^.r). (2) 

h,k 

In the Abrikosov solution of the Ginzburg-Landau (GL) 



equations (as applied to a square lattice) Jig] , the Fhk 
are given by: 

F hk oc _(_i)(^ 2 +^) . ex p ( - |(ft 2 + fc 2 )) ; (3) 

this rapidly falls off with q (see Table |). 

The Abrikosov solution is only valid near B C 2 . In high- 
k superconductors, with the field not close to B C 2, the 
London expression [|l9| is appropriate instead. This gives 
F hk oc 1/(1 + q\ k \ 2 ). Note that unlike the Abrikosov 
solution, all the Fhk are positive. Table Q shows that 
the Fourier components fall off much less rapidly with 
q. However, strontium ruthenate has a value of the 
Ginzburg Landau parameter k = A/£ ~ 2.0 for the field 
along the c axis, which means that the London approach 
is not realistic except at very low inductions. Therefore, 
to see what conventional GL theory predicts for this ma- 
terial at lower fields, one must use the Brandt numerical 
solution of the GL equations [ p0| . Typical results are 
given in Table [l[ 

Next, we consider the Agterberg TCGL solution 
which is equivalent to the Abrikosov one, except that 
there are two complex order parameters instead of only 
one. In the mixed state with B parallel to c both com- 
ponents are automatically present because of mixed gra- 
dient terms in the free energy functional (|J§]. Typical 
values from this theory for Fhk, relative to F\o and the 
resulting SANS intensities are given in Table H It may 
be seen that the two-component theory gives intensities 
that fall off much less rapidly with q than those given by 
the one-component Abrikosov solution. 

Under the conditions of our experiments, where the 
field is not close to £? C 2, it may be argued that the 
Abrikosov approximation used by Agterberg is not ap- 
propriate. However, recently Heeb and Agterberg have 
solved numerically the GL equations at all fields for the 
TCGL case pl|| . We also give in Table | a list of Fourier 
components from these calculations, using values of pa- 
rameters that appear to describe our results quite well. 

The corresponding vortex structures in real space are 
shown in Figures Note that there is a minimum field 
point in the two-component theory (for the conditions of 
our experiment) which lies between the positions of the 
flux line cores, not in the centre of the square. We give 
results of this theory for two values of the parameter v 
(— 1 < v < 1) |9|] which describes the degree of fourfold 
anisotropy of superconducting electrons (y = corre- 
sponds to a cylindrical Fermi surface). We note that the 
results do not change greatly with v. Hence the quali- 
tative difference between Figures and || is due to the 
difference between TCGL and GL theories rather than 
effects of fourfold crystal anisotropy. It may be that v 
is quite small since \v\ > 0.0114 is sufficient to stabilise 
a square FLL and align it to the crystal lattice with an 
orientation determined by the sign oiv ||] . 

We now turn to measurements of the FLL structure. 
Single crystal Sr2Ru04 was prepared by the floating zone 
technique with excess Ru02 as a flux |23]. Six plates of 



total mass 556 mg were cleaved from the as-grown crys- 
tal and annealed for 72 hours in air at 1420° C to remove 
defects and increase T c , which was 1.39K with a width 
(10-90%) of ^50 mK. With the field applied parallel to 
the c-axis at 100 mK, the value of B C 2 was 58mT. For the 
small angle neutron scattering (SANS) measurements, 
the samples were mounted with conducting silver paint 
as an aligned mosaic with their c-axes perpendicular to 
a copper plate, which was mounted on the mixing cham- 
ber of a dilution refrigerator. This was placed between 
the poles of an electromagnet, which had holes parallel to 
the field for transmission of neutrons. The magnetic field 
was parallel to the c-axes of the crystals within 2°, and 
the FLL was observed using long-wavelength neutrons 
on instrument D22 at the Institut Laue Langevin. Typi- 
cal wavelengths employed were 14.6 A, with a wavelength 
spread (FWHM) of 12%; the neutron beam was incident 
nearly parallel to the applied field, and the transmitted 
neutrons were registered at a 128 x 128 pixel multidetec- 
tor (pixel size 7.5 x 7.5 mm 2 ) placed 17.71m beyond the 
sample. Typical results are shown in Fig ||. In addition 
to the strongest {10} reflections, the {11} reflections are 
strong, and higher orders are present. The intensity of 
the strongest diffraction spot is <10 -3 of the incident 
beam intensity, so these higher order reflections are not 
due to multiple scattering. Their intensities are recorded 
in Table |: it will be noted that they are much larger 
than those given by the Abrikosov structure. 

To reconstruct the B(r) of the FLL corresponding to 
these results, we require the sign of Fhk relative to Fio 
(the FLL is centrosymmetric, so all the Fhk are real). 
The most important component after Fio is Fu. If it 
has the same sign as Fio, then the {11} components add 
in phase at the flux line cores to give a field peak that 
is sharper than the field minimum. Measurements of the 
field distribution in strontium ruthenate by /iSR [ p3|]24| 
show that this is the case. This sign for Fu is not sur- 
prising, since all models in Table Q give it as positive. For 
the small contributions of F20 and F21, we may assume 
the same signs as given by the Agterberg and Abrikosov 
solutions: taking the London sign makes a large differ- 
ence to B(r), and also can be ruled out by fiSR results . 
The reconstruction of B(r) is shown in Fig. 4. Note that 
it is completely different from the Abrikosov or Brandt 
solutions to the GL equations, and in good qualitative 
agreement with the TCGL predictions. 

The results we have given so far correspond to low 
temperature and a particular magnetic field value. In 
Table || we present the values of the form factors Fio and 
Fu for a range of fields at lOOmK. Also, in Fig 3 we plot 
versus temperature the ratio of the Fourier components 
for the strongest two reflections Fn/Fio at 10, 20 and 
30mT. Remarkably, this ratio varies little with field and 
temperature and does not tend to the Abrikosov value 
as T^T C . Non-local effects JHJ and deviations from GL 
theory in ultra pure superconductors (2^| should both die 
away at high temperatures. Therefore, these effects are 
not expected to be the cause of the flux line shapes we 
report, although they may affect the details of B(i) at 



low temperatures. 

In conclusion, the strength of the higher order reflec- 
tions from the FLL in strontium ruthenate and their tem- 
perature dependence certainly show that a standard one 
component Ginzburg Landau model is insufficient to ex- 
plain the observed diffraction pattern. However, our re- 
sults are in good qualitative but not perfect agreement 
with a two component Ginzburg Landau theory. Uncon- 
ventional flux line shapes in this material are strong evi- 
dence for unconventional superconductivity in Sr2Ru04. 
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ported by the U.K. E.P.S.R.C., CREST of Japan Science 
and Technology Corporation, and the neutron scattering 
was carried out at the Institut Laue-Langevin, Grenoble. 
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FIG. 1. Contour plot of the magnetic field in a square 
flux line lattice as given by Brandt's numerical solution of 
the Ginzburg-Landau equations for the particular case of 
k — 2.0, B — 20raT and B C 2 = 58raT. Contour lines are 
equally spaced. This result is very similar to that given by 
the Abrikosov solution valid near H C 2 1181. 
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FIG. 2. Contour plot of the magnetic field in a square flux 
line lattice as given by Agterberg's solution of his TCGL equa- 
tions, valid near B C 2, in the case of a cylindrical Fermi surface 
(i/ = 0.0). 
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FIG. 3. As for Figure u\ but with a fourfold distortion to 
the Fermi surface (y — 0.2). 
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FIG. 4. Heeb and Agterberg's numerical solution to the 
TCGL equations with the parameters k = 1.6, v — 0.05 
chosen to give a good fit to our data, with an applied field of 
B 20mT and £ c2 (100mK)=58mT. 
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FIG. 5. Contour plot of FLL diffraction pattern. A field of 
20mT was applied parallel to c above T c ; the weak diffracted 
beams due to the flux lattice were extracted from background 
scattering by subtracting data taken above T c . The axes are 
pixel numbers and the central region of the detector has been 
masked. The data shown is a sum of 5 patterns obtained by 
rocking the FLL up ±0.3° about the [110] axis. 
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FIG. 6. Temperature and field dependence of the ratio 
Fu/Fio. The {1,1} intensity used is the direct average of 
all four spots, and the {1,0} intensity is a weighted average 
of the top and side spots which allows for the different rock- 
ing- curve width in the vertical and horizontal directions Jl^] 
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FIG. 7. Contour plot of the magnetic field in the mixed 
state of SRO as reconstructed from the data represented in 
Figure ^, using the signs of the Fourier components as dis- 
cussed in the text. 



TABLE I. Calculated and experimental flux lattice Fourier 
components and intensities for B=20mT, B C 2=58mT . 



/i, k of diffraction peak: 


{1,1} 


{2,0} 


{2,1} 


qhk/ qio 


1.414 


2.0 


2.236 


Fhk/Fio (London A = 152nra) : 


0.53 


0.27 


0.22 


Fhk/Fio (Abrikosov): 


0.2079 


-0.00898 


0.00187 


F hk /F 10 (Agterberg v = 0.0): 


0.5657 


-0.1051 


0.0353 


F hk /F 10 (Agterberg v = 0.2): 


0.7711 


-0.0667 


0.0503 


Fhk/Fio (Heeb and Agterberg 


0.484 


-0.019 


0.046 


z/ = 0.05, k = 1.6): 








Ihk/ho (London): 


0.199 


0.0365 


0.0216 


Ihk/ho (Abrikosov): 


0.0306 


0.00004 


0.000002 


Ihk/ho (GL (Brandt)): 


0.0783 


0.00091 


0.000298 


Ihk/ho (Agterberg v = 0.0): 


0.2263 


0.00552 


0.00056 


Ihk/ho (Heeb and Agterberg 


0.166 


0.00095 


0.00018 


v = 0.05, k = 1.6): 








Ihk/ ho 


0.197(2) 


0.011(3) 


0.007(1) 


(Expt. at 20 mT lOOmK): 









TABLE II. Experimental FLL Fourier components at 
lOOmK. Errors represent statistical errors; FLL disorder or 
static Debye Waller factors may reduce Fhk for high q. 



Magnetic Field (mT) 



F10 (mT) 



F11 (mT) 



100 
150 
200 
250 
300 
350 
400 



1.14(3) 
1.02(9) 
0.88(1) 
0.71(1) 
0.58(1) 
0.45(1) 
0.30(2) 



0.66(3) 
0.54(1) 
0.47(1) 
0.36(1) 
0.29(1) 
0.24(1) 
0.14(1) 



